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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE 3725 


AERODYNAMIC INTERFERENCE OF SLENDER WING-TAIL 
COMBINATIONS 
By Alvin H. Sacks 


SUMMARY 


Mathematical expressions are derived for the interference forces 
and moments acting on the tails of slender plane and cruciform wing-tail 
combinations of general plan form in steady straight flight at combined 
angles of attack and sideslip. The derivations are made within the limi- 
tations of slender-body theory under the assumption that the vortex sheet 
leaves the wing as a flat sheet and becomes fully rolled up ahead of the 
tail. The derived expressions are used to calculate the steady lifts, 
side forces, pitching moments, and rolling moments of a number of wing- 
tail combinations. The effects of changes in tail height, tail length, 
ratio of tail span to wing span, tail incidence, and tail thickness are 
calculated. The resulting curves , and particularly their nonlinearities, 
are discussed at some length in connection with static stability. In 
general, the most dramatic effects are noted when the vortices shed from 
the wing strike the tips of the tail trailing edge. 

An expression is developed for the lift of a plane wing-tail combina- 
tion which is pitching and plunging, and the associated stability deriva- 
tives are calculated as functions of the angle of attack. Discontinuities 
in lift-curve slope C L and the stability derivative C L& are noted 
for plane wing-tail combinations with high tails if the span of the tail 
is slightly greater than the span of the vortices shed by the wing. 

Photographs of the wake of the wing in the presence of the tail, as 
observed in a water tank, are presented for a plane triangular wing- tail 
combination with a high tail. The measured variation of the lateral 
spacing of the wing vortices with distance in the presence of the tail is 
presented an d discussed in connection with the assumptions of the analysis. 
The tail was found to cause an appreciable inboard shift of the wing vor- 
tices for the case investigated. 

INTRODUCTION 


It is now generally recognized that for wing-tail combinations 
involving low-aspect-ratio wings and moderate tail lengths it is not 
permissible to neglect the rolling up of the vortex sheet in calculating 
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wing-tail interference effects. Since such cases are often encountered 
in connection with high-speed airplane and missile designs, some effort 
has been devoted in recent years to the calculation of wing-tail inter- 
ference, accounting in some manner for the rolling up of the vortex sheets 

Except for some numerical work in tracing the rolling-up process 
itself (e.g., refs. 1 , 2 , and 3 )> the distortion of the vortex sheet is 
usually accounted for simply by assuming the sheet to be fully rolled up, 
and the emphasis has generally been on obtaining expressions for the 
forces on the tail in terms of the positions of the vortices shed by the 
wing (e.g., refs. 4, 5 , 6 , and 7)* In 1948, however, Graham (ref. 8 ) 
actually calculated the variation of tail lift with angle of attack for 
some planar wing-tail combinations. This was done by expressing the 
vortex positions as a function of the angle of attack under the two alter 
native assumptions of a flat vortex sheet and a fully rolled-up vortex 
sheet shed by the wing. Graham's results, which for the slender tail 
case were only approximate , showed some interesting departures from the 
usually assumed linear variations. In 1952 Morikawa (ref. 9 ) investi- 
gated the "maximum" wing-body-tail interference (which occurs when the 
rolled-up vortices lie in the plane of the tail) by restricting the 
analysis to angles of attack near zero and at the same time assuming that 
the vortex sheets are fully rolled up. This avoids any consideration of 
the nonlinear variation of tail lift with angle of attack. 

In the present paper, the emphasis is placed on calculating the 
variations of total forces and moments with angles of attack and sideslip 
for a number of slender plane and cruciform wing-tail combinations and 
for some airplane-type arrangements of a plane wing and a horizontal and 
vertical tail. Significant nonlinearities are found, and these will be 
discussed in some detail with regard to their effects on the static 
stability of the various combinations. The lift of a plane wing-tail 
combination which is pitching and plunging will also be determined and 
the variation of pertinent stability derivatives with angle of attack 
will be calculated. 

The primary assumption in the present analysis will be that the 
vortex sheet leaves the slender wing as a flat sheet and becomes fully 
rolled up ahead of the tail. It will also be assumed, as is customary, 
that the tail does not influence the positions of the vortices shed by 
the wing. The validity of both of these assumptions has been investi- 
gated experimentally by means of a water tank, and these results will 
also be presented and discussed. 


SYMBOLS 


A aspect ratio 

b local semithickness of horizontal tail 
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maximum chord of wing 

distance from wing apex to pivot point 

maximum chord of tail 

tail length, Z-c 

distance behind wing trailing edge at which vortex sheet is essen- 
tially rolled up 

height of trailing edge of horizontal tail above wing chord plane 
incidence of horizontal tail relative to wing chord plane, radians 
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force in the z direction (approximately lift) 

rolling moment about x axis 

over-all length of wing-tail combination 

number of external (free) vortices 

pitching moment about pivot point x = c x 

yawing moment about pivot point x = c x 

angular rolling velocity about the x axis, radians/sec 

angular pitching velocity about the pivot point x = c x , 
radians /sec 

V + iW 

angular yawing velocity about the pivot point x = c x , radians/sec 

radius of transformed circle corresponding to wing or tail cross 
section 

local semi span of wing or tail 
airplane cross-sectional area 

plan-form area of wing (area of one wing of cruciform) 
maximum semi span of wing (at x = c) 
maximum semispan of tail (at x = l) 
time, sec 

maximum span of vertical tail panel 

local span of upper vertical tail panel 

local span of lower vertical tail panel 

component of flight velocity along the negative x axis 

component of flight velocity along the positive y axis 
(V 0 = U 0 3 if P = 0) 

v 0 - r ( x - c x ) 

component of flight velocity along the positive z axis 
(W Q = -U Q a if P = 0) 
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' w 

0 * 

1 

ST 

1 

0 

H 

• X TE+ 

distance from wing apex to position immediately behind wing 
trailing edge 

xyz 

Cartesian coordinates fixed in the tail as illustrated in 
figure l(a) 

Xl 

distance behind the wing trailing edge 

x' 

distance behind apex of the tail, x x - d + c' 

Y 

force in the y direction (side force) 

y i> z i 

y and z coordinates of starboard rolled-up vortex 

a 

angle of attack, radians 

a cr 

angle of attack at which the rolled-up wing vortices intersect 
the line containing the trailing edge of the horizontal tail 

a, 

time rate of change of angle of attack, radians/sec 

P 

angle of sideslip, radians 

r 

strength of one rolled-up vortex shed from the wing 

r k 

circulation of kth external (free) vortex, positive counter- 
clockwise 


y + iz 

^k 

position of kth external (free) vortex, y k + iz k 

i 

*1 

s o 

P 

fluid mass density 

a 

complex coordinate in transformed circle plane 

°k 

position of kth external vortex in transformed circle plane 

a *r 

position of kth external vortex relative to its image in the 

r 2 

transformed circle, cr v - — — 

°k 



b 

s 

horizontal-tail thickness ratio (constant for conical tails treated 
herein) 
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SUBSCRIPTS 


I due to vortex interference 

l tail trailing edge 

T tail (when used on forces and moments indicates isolated tail) 

TE wing trailing edge 

w wing 


SPECIAL NOTATIONS 



contour integral taken once round the cross section in the positive 
(counterclockwise) sense 


R real part 


I imaginary part 

( ) complex conjugate of ( ) 


ANALYSIS 


The present report is concerned with the calculation of the total 
aerodynamic forces and moments exerted on some plane and cruciform slender 
wing-tail combinations in steady and maneuvering flight. The calculations 
will be made within the limitations of slender -body theory and will employ 
the techniques developed in references 10 and 11. Inasmuch as it was 
shown in the former reference that the forces and moments are linear in 
the potential (although not in the motions), it is permissible to calculate 
the forces and moments due to vortex interference alone and add them to 
those of the isolated wing and tail. Consequently, the major portion of 
this report will be concerned with the calculation of the forces and 
moments due to wing-tail interference. 


It was shown in reference 11 that the components of force and moment 
due to vortex interference can be expressed as 


tj - iLj = -ipU Q 
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( 1 ) 


MCA TN 3725 


7 



( 2 ) 



(3) 


where R denotes the real part, p is the rate of roll about the x 
axis, F' is the complex potential due to the wing vortices and their 
images in the tail, and a^ r represents the (complex) distance between 
the kth shed vortex and its image in the transformed circle plane. 

Evidently, the two essential quantities to be determined before equa- 
tions (l) to (3) can be applied are the additional complex potential F' 


their images . 

Although the analyses of references 10 and 11 employed a coordinate 
system whose x axis passes through the airplane nose , it can be shown 
that the results of those analyses are unaffected by a normal translation 
of the x axis provided that the cross-sectional area satisfies the condi- 
tion S = (ds/dx) = 0 at the airplane nose. This condition is satisfied 
by 1 the wings of the wing-tail combinations to be treated in the present 
analysis, and it is convenient for our purpose here to use a coordinate 
system fixed in the tail as illustrated in figure l(a) . 

It can be seen from equation (l) that for steady straight flight 
[P = (S/St) = 0] the interference lift and side force are independent of 
the plan forms of the wing and tail. Consequently, the calculated curves 
of steady lift and side force presented in this report are for slender 
wings and tails of arbitrary plan form except that their trailing edges 


and the sum 


k=i 



representing the impulse of the shed vortices and 
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must lie in planes normal to the x axis. On the other hand, equa- 
tions (l) , (2) , and ( 3 ) show that the lift and side force in unsteady- 
flight and the moments in steady flight require integrations over speci- 
fied plan forms. The calculations of steady moments and of unsteady 
lift presented in this report have been carried out for triangular plan 
forms only. 


STEADY STRAIGHT FLIGHT 


In this portion of the analysis we shall specialize equations (l) 
and ( 2 ) to steady straight flight by setting p = (8/dt) = 0 . The last 
two terms of each of these equations are thus eliminated. 


Plane Wing and Cruciform Tail 


For most of the wing-tail combinations to be treated here, the wing 
is taken as a planar, slender, pointed thin wing having its maximum span 
at the trailing edge which is assumed straight. In all cases, the tail 
is considered rigidly attached to the wing. It is also assumed that the 
vortex sheet leaves the wing as a flat sheet at the trailing edge and 
becomes fully rolled up into two discrete line vortices somewhere ahead 
of the tail. 


It is further assumed that the tail does not influence the vortex 
positions. The validity of this assumption has been investigated experi- 
mentally and will be discussed in a later section of this report. For 



such cases, then it can be seen from sketch (a) that the vertical posi- 
tion z x of both vortices relative to the horizontal tail is given by 


zi = (a w - e)xx - hip - (d - Xi)i T (4) 


where the angle e is easily computed by calculating the velocity of a 
two-dimensional vortex pair placed at the centroids of vorticity of the 
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wing. Since this procedure places the vortices (jt/2)s 0 apart, due to 
the elliptic circulation distribution (which is unaffected by sideslip) 
the angle e is found to be equal to (2/ J t 2 )a v . It should be mentioned 
that the positions calculated in this way are in good agreement with 
those observed experimentally behind a wing alone, as pointed out in 
reference 12, even though the details of the rolling up have been ignored. 
The lateral positions (in body axes) are evidently 

yi = f s o - 3 x i > y 2 = - f s o - P X 1 (5) 

where the subscripts 1 and 2 on y refer to the starboard and port vor- 
tices, respectively. Thus the positions of the two rolled-up vortices 
at a distance x x behind the wing are completely defined by 


£1 


£ 2 


y x + iz x = ^ s Q - i(h T + diij) - 
y 2 + iz-L = - s Q - i(h T + di T ) - 


x i{p -i[aw(l - Jr) + } 

xxjp - i[av(l - |f) + It] } 


(6) 


The magnitude of the strength T of each of the rolled-up vortices 
is obtained by equating the impulse of the vortex pair to the lift of 
the wing , that is , 

pU o r 75 s Q = Lv = JtpU Q Ov/Sq 2 


or 

f = SUqO^-Sq 


(7) 


It is noted that for the fully rolled-up vortex pair m = 2 and 
r 2 = ~T X = -r (see sketch (b)). Finally, in order to calculate the 



Sketch (b) 
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m 

sum ^ the transformation of the tail cross section to a circle 


k=i 


(leaving the flow field at infinity unchanged) is required. Such a 
transformation for the class of cross sections shown in sketch (b) is 
given in reference 10 as 


t 2 - 


= [! 


(h - f) + a - 


(h + ff 
l6a 


( 8 ) 


and 


h + f 

4 


where 


h = o/s 2 + t x 2 ; f = J s 2 + t 2 2 


so that the inverse transformation relating the vortex positions in the 
two planes is 



/ r / \ r q 2 I 

r q2 

/ [ i (Ts-t t! 2 - Js* + t a «j - 2 7^ k 2 - S 2 1 + 

k/s 2 +t 1 2 +Vs 2 +t 2 2 


( 9 ) 


In writing the transformation in this form, we require that square roots 
of complex quantities be evaluated in the manner of the appendix. The 
proper signs are thus automatically taken care of. 

m 

The Siam ^ r k a k r can now be determined from equations ( 7 ) and ( 9 ) 
k=i 

if it is recalled that 

% s ^ < 10 > 

With the expression for this sum, then, all the force and moment compo- 
nents except the rolling moment can be calculated directly from equa- 
tions (l) and (2) for the wing-tail combinations discussed above. It 
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m 

will be noted that the above expressions for the sum ^ ^k g k r have been 

k=i 

developed only for the fully rolled-up vortices. However, between the 

m 

wing and tail, since there is no body, the impulse pU Q r k cr kr is inde- 

k=i 

pendent of the wake shape , being always equal to the lift of the wing. 
Therefore, it has only been assumed that the vortices become fully rolled 
up somewhere ahead of the tail. 


The following paragraphs will be concerned with the actual calcula- 
tions relating to the longitudinal , directional, and lateral stability of 
several wing-tail combinations falling into the category discussed in the 

foregoing analysis. The determination of the integral J^F '&((£) will be 


carried out in the discussion of lateral stability since it arises only 
in the calculation of the rolling moments . 


Longitudinal stability . - Since the vertical tail has no effect on the 
purely longitudinal 0 = 0) aerodynamic characteristics of the wing-tail 
combination, we can for this discussion set t^ = t^ = 0, thereby simpli- 
fying the calculations. In fact, since equation (8) then reduces to the 
Joukowsky transformation, the addition of elliptical thickness to the 
horizontal tail offers no difficulty and will be incorporated. Thus the 
expression for to be used for calculating the longitudinal charac- 

teristics will be (see sketch (c)) 




Sketch (c) 
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°k = | [c k + JC k 2 - (s 2 - b 2 )j 


(ID 


and 


r o = 


s + b 


The thickness of the tail will be kept small in the calculations, in view 
of the assumption that the tail does not affect the vortex positions. 

The wing-tail combinations treated in this section will have the general 
appearance shown in figure l(b) . The wing is shown as a triangular flat 
plate and the horizontal tail is shown as a thin elliptic cone since these 
assumptions will be made in the calculation of the pitching moments . 
However, they will not be made for the calculation of the steady lift 
since equation (l) shows that the interference lift and side force in 
steady straight flight are independent of plan form. 

If it is noted that, due to symmetry (see sketch (c)) cr 2 = - 0 1} the 


sum 


£r k a kr over the tail can be written directly as 


k=i 


C- 

l 

k=i 


r k O kr , = r i a i r + r 2 a 2 r 


= r (°i r + 5 i r ) = 2rR ( ff i r ) 


( 12 ) 


But from equations (10) and (ll) we have 


- £ 


- 1 {?! + Ji ? - - * 2 > - (H-£)[si - AY - <= 2 - b2 >] } (13) 

so equation (12) becomes 


= tr + A/ - U 2 - b 2 ) - 


k=i 


(s 2 - b 2 )l j 
(14) 
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Returning now to equation (l) we have for steady straight flight 



and, by virtue of symmetry for 3 = 0, the interference side force Yj 
vanishes. (This also follows from eq. (l4) which shows that the above 
stuns are real . ) The interference lift Lj is given by the negative of 
the imaginary part of equation ( 15 ) • 



which can now be evaluated from equations (6), ( 7 ), and (l4) . That is, 
from equation (6), = (n/4)s 0 for all x and we note that imme- 

diately behind the wing s = b =0. Thus it follows from equations ( 7 ), 
(l4) , and ( 16 ) that 


L I 


4 

1- (b/s) z 


PU 0 o-wSo 




(17) 


where s x is the semispan of the tail at its trailing edge. The quan- 
tity C, x gives the position of the starboard vortex at x = l and is 

evaluated from equation (6) by setting 3=0 and x x = d. The total lift 
of the wing-tail combination is obtained by adding to equation ( 17 ) the 
lifts of the wing and tail alone as given by slender -body theory (see 
eq. ( 7 )). The resulting expression for the lift coefficient, based on 
the area of the wing, is 


Cl 


= | Aw(c 


CL\j + <Xrp 



2AvOv 

1- (b/s)j 





n 

4 


(18) 


and it can be seen that, to the order of the present analysis, tail 
thickness affects the lift only if the tail has a blunt trailing edge. 
The procedure for calculating the real part of the square root without 
ambiguity has been discussed in Appendix B of reference 11. 
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By making use of the same substitutions in equation (2) as were used 
in equation (l) for the interference side force and lift, one finds simi- 
larly that the interference yawing moment vanishes (for 3=0) and that, 
after integration by parts, the interference pitching moment about the 
pivot point x = c x can be expressed as 



provided that b/s is constant over the tail. 1 In this expression c ! 
is the chord of the tail and x ! is distance along the x axis measured 
from the apex of the tail. Now let us choose the pivot point to be 
c 1 = 2/3 c and the plan forms of wing and tail to be triangular. Then 
the triangular wing contributes no pitching moment, and if we add the 
known pitching moment of the conical tail alone as given by 

m t = -L T (a + £-f.) = - Ttpu/oTS! 2 (a + |-^) (20) 


the resulting expression for the total pitching-moment coefficient, based 
on the wing area s Q c and the wing chord c is 


^m - 




2Ava w 
1 - (b/s) 


(t + i) R 


t 2 

so 2 



- + - 




( 21 ) 


The integration indicated in equation (21) is most easily carried out if 
the real part is taken after the integration. With this in mind, then, 
we observe that the variable x f is related to x± by (see fig. 1(b)) 

1 The corresponding expression for arbitrary chordwise thickness 
distribution (b/s f const.) can just as easily be obtained in the same 
way, but the calculated results will depend on the particular shape chosen 
and the integrals will in general have to be evaluated by numerical or 
graphical methods for each case . 
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Xl = x 1 + d - c f (22) 

and that for the conical tail s/s Q = (A T /4) (x f /s 0 ) and b/s is constant. 
Hence } by means of equation (6), the complex square root in the integrand 
can be expressed as 



and the integral in equation (21) can be expressed in the form 


J S ° H Jx dt, = IE (2D! + C)\/x + ~ ln(jx + + 


2n/d/_ 


1 _ 

4d 


R (2D| + C) R -Jx ~ I (2D| + C) I n/T + 






( 25 ) 


where | = ( x'/s Q ). Now if any of the terms in this expression has a 
discontinuity within the interval of integration the integral must be 
divided at that point. Such discontinuities will in fact occur if the 
rolled-up vortex passes through the horizontal tail. This occurrence 
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will be discussed later from the physical point of view, but for the 
moment we need only to recognize that the mathematical expression of this 
situation is (for a triangular tail) 


H - 0 for JL < (26) 

a T b o s o 

where z ± is given by the imaginary part of equation (6) . The variations 
of Rjx and I-Jx for a typical case in which equation ( 26 ) is satisfied 
are shown in sketch (d) and it can be seen that at g = g (the value of 
x'/ s o at which z x is zero) the real part of the square root vanishes 
and the imaginary part has a discontinuity. 



Thus, making use of the continuity of the real part of the square 
root at |=|, and writing out the real and imaginary parts of the 


logarithm, we can finally write equation (25) in the form 


C | 

J S ° U -JX di = J Rn/X d£ + J S ° Rn/x 


d£ 





Zn 


n/F + i«/d + 



- I(2DS + C) 





where G is identical with the quantity in the brackets of equation (25) and the subscripts 
£ 0 - and i Q + refer, respectively, to values immediately before and after the discontinuity. The value 
of n (an integer or zero) is determined by the fact that the resulting expression for the integral 
must be a continuous function of since B, C, and D are continuous functions of a^. 


Calculations have been carried out by means of equation (l8) and equations (21) through (26a) to 
determine the total lift and pitching-moment coefficients of a number of wing-tail combinations of the 
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type shown in figure l(b) . It can he seen from equations (l8) and (21) 
(recalling eq. (o)) that the theoretical lift and pitching moment are 
nonlinear functions of the angle of attack c^. In view of this fact, 
it seems appropriate to postpone temporarily the calculation of stability 
derivatives and to calculate instead the actual variations of forces and 
moments with angle of attack. For the present calculations, the aspect 
ratio of both wing and tail was chosen as 2 in all cases, and the effects 
of horizontal tail height hip, tail-span-to-wing-span ratio s 1 /s 0 , tail 
incidence ip, horizontal-tail thickness ratio b/s, and tail length d 
have been investigated. The calculated results are presented in figures 2 
to 6 and the lift coefficients are divided by (n/2)A w since equation (l8) 
shows the lift coefficient to be linear 2 with respect to the wing aspect 
ratio. A w . The pitching moments are nonlinear in Ay since equation (21) 
contains products of Ay with d/c and c'/c which are themselves propor- 
tional to A w since all quantities are specified in terms of the wing 
semi span s Q . 

The effect of horizontal -tail height on the lift and pitching moments 
is illustrated in figure 2 for zero tail incidence and zero tail thickness. 
It can be seen from equation (17) that for b/s = 0, the interferenc e lift 
becomes equal and opposite to the lift of the wing if - s x 2 = 0. 

This condition is satisfied if £ is real and less than s x , that is, 

if the rolled-up vortices from the wing just intersect the trailing edge 
of the tail. In figure 2, the tail heights have been chosen so that this 
situation occurs at simple values of a w , and it can be seen that, in 
each case, at that angle of attack the total lift is equal to the lift of 
the tail alone which is equal to the lift of the wing since the tail span 
is equal to the wing span. It would appear from figure 2 that unless the 
horizontal tail is placed either low (hp < 0) or very high, the lift and 
pitching -moment curves are apt to be highly nonlinear, even in the rela- • 
tively low positive angles of -attack range. Severe static instability 
precedes the "critical" angle of attack (the angle of attack at which the 
wing vortices intersect the tail trailing edge) and abrupt changes in 
both lift-curve slope and static stability are observed at the critical 
angle . 

While the sharp breaks in the lift and pitching -moment curves indi- 
cated by the theory will no doubt be smoothed out somewhat in actual 
flight through a real fluid, it might be well to consider the signifi- 
cance of the critical angle of attack described above. First, the 
existence of such a well-defined critical angle is predicated on the 
assumption that the vortices shed from the wing lie in a horizontal line 
when viewed in planes x = const., that is, that the vortex sheet either 
remains flat or is fully rolled up ahead of the tail. In such a case, 
the tail experiences a maximum download when the line containing its 

It must be noted that the equations also show the lift and pitching 
moment to be nonlinear with respect to all of the parameters investigated 
in the calculations. The calculated effects of the separate parameters 
are therefore not additive. 
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trailing edge intersects the vortices, as the tail is then in a maximum 
downwash field. Such a critical angle might be thought of as a local 
stall angle for the airplane and will be predicted by the theory for all 
slender plane wing-tail combinations whose wing and tail do not lie in 
the same plane. It might be said that if the wing and tail lie in the 
same plane, the critical angle is zero, but since the downwash is then 
also zero, the angle has little significance. 


The critical angle of attack can be calculated directly from equa- 
tion (6) by setting z x (i.e., the imaginary part of or £ 2 ) to zero 
at x x = d. Thus one finds 


Ocr - 


hj/sp 

(d/s 0 )[l- ( 2 /fl 2 )] 


(27) 


for the combinations treated in this section. Evidently, then, <xc r 
depends on the ratio of tail height to tail length and will be positive 
if the tail is above the plane of the wing and negative if the tail lies 
below the plane of the wing. For the former case (high tail), the tail 
will experience a maximum download, causing a peak positive pitching 
moment, while for the latter case (low tail), the tail will experience a 
maximum upload, causing a peak negative pitching moment. These two types 
of cases are clearly seen in figure 2. The critical angle of attack can 
of course be made larger than the stall angle of the wing itself by 
placing the tail very high. In figure 2 for hj/sp = I.9137, for example 
the lift and pitching-moment curves remain fairly linear up to an angle 
of attack of about 10°. It should be mentioned with regard to the above 
discussion that, if the vortex sheet is only partially rolled up at the 
tail location, the vortices do not lie in a horizontal line there, the 
critical angle of attack loses definition, and the lift and moment curves 
although nonlinear, will not exhibit the sharp breaks discussed above. 


In figure 3> the ratio of tail span to wing span is varied for a 
given tail height. It can be seen that the lift and pitching-moment 
curves exhibit sharp peaks when the ratio of tail span to wing span 
becomes equal to or slightly greater than jt /b (the ratio of vortex span 
to wing span). That is, once again, the curves have sudden changes in 
slope if the wing vortices strike the tail trailing edge near the tips. 
Thus, it might be said that there is a critical ratio of tail span to 
wing span as well as a critical angle of attack. However, it should be 
pointed out that the specific value of tc/ 4 is a consequence of several 
assumptions; namely, (l) that the vortices are fully rolled up ahead of 
the tail, (2) that the tail has no influence on the vortex positions, 
and (3) that the trailing-edge cross section of the tail is f ull y effec- 
tive in destroying the downwash field of the wing, even though the chord 
of the tail vanishes at the tips. This last assumption is associated 
with the slender-body assumption that the flow in planes x = const, 
is essentially two-dimensional. 


The sign of the pitching moment at the critical angle of attack can 
be determined from the ratio of tail span to wing span, provided this 
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ratio is equal to or greater than the critical ratio. It has already- 
been pointed out that the lift of the plane wing-tail combination at the 
critical angle of attack is equal to the lift of the tail alone at its 
geometric angle of attack, if the tail span is at least as great as the 
vortex span. Therefore, the lift on the tail of the combination is equal 
to the difference between the lifts of the isolated tail and wing at 
their geometric angles of attack. But, according to slender-wing theory, 
these lifts depend only on the maximum spans of the wing and tail. Hence, 
it can be concluded that if the span of the tail is less than the span of 
the wing, the lift on the tail at the critical angle of attack is negative 
and the pitching moment of the combination is therefore positive. Since 
it was stipulated above that the ratio of tail span to wing span was 

greater than the critical value, the conclusion is that for ^ <1, 

the pitching moment at the critical angle of attack is positive . 8n the 
other hand, if the tail span is greater than the wing span, the lift on 
the tail is positive at a = acr> so the pitching moment of the combina- 
tion is negative. Finally, if sx/so = 1> the lift on the tail, and 
hence the pitching moment of the combination, vanishes at the critical 
angle of attack. 


The above conclusions may be summarized schematically as shown in 
sketch (e) . It can be seen that in the first case there are three trim 



Sketch (e) 


points, two stable, and one unstable; in the second case there is only 
one trim point, which is stable; and in the third case there are two trim 
points, one stable and one neutrally stable. In the above discussion, 
the tail incidence has been assumed to be zero. It has also been assumed 
that the pitching moment of the wing alone is zero, which is the case for 
all the combinations treated in this paper. 


The effect of changes in tail length on the lift and pitching-moment 
curves is shown in figure k for two tail heights with the tail at zero 
incidence and the tail span equal to the wing span. For zero tail height 
(fig. 4(a) ) , the predominant effect of increasing the tail length is 
apparently the larger pitching moments due to the increased lever arm. 
However, the lift on the tail at a given angle of attack is also seen to 
increase with the tail length. This is due to the greater vertical dis- 
tance between the tail and the wing vortices, which results in a reduction 
of wing-tail interference. 
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For the case of a given nonzero tail height (fig. 4 (b), h T /s o = 0 . 9568 ) , 
the most obvious effect of a change in tail length is a corresponding 
change in the critical angle of attack. It can be seen that reducing the 
tail length ratio d/s Q from 6 to 4 increases from 0.20 to 0.30, 

as predicted by equation (27). Now it is interesting to compare fig- 
ure 4 (b) with figure 2 in which a cr was changed by a change in tail 
height with a fixed tail length. Let us say, for example, that we start 
with a wing- tail combination having d/s Q = 6 and a C v. = 0.2, and we wish 
to shift oc r to 0 . 15 . A comparison of figures 4 (b) and 2 shows that, 
if this shift is accomplished by increasing the tail length ratio to 8, 
the resulting pitching moments are about 50 percent larger than if it is 
accomplished by lowering the horizontal tail. 

In figure 5 , the tail incidence is varied at two different tail 
heights for a fixed tail length (d/s 0 = 6) and a fixed span ratio 
(si/s 0 = 1 )- It can be seen that for hj = 0 (fig. 5 (a)), the lift and 
pitching-moment curves are shifted linearly up and down with tail inci- 
dence. Since the curves are slightly nonlinear in the angle of attack, 
however , the result is a slightly nonlinear variation of trim lift coef- 
ficient with control setting (for an all -movable tail) . Now if the tail 
does not lie in the plane of the wing, (h T /s 0 = O.9568, fig. 5(b)), the 
nonlinearities in the lift and pitching -moment curves, and hence in 
control effectiveness, become more pronounced. It can be seen in fig- 
ure 5(b) that the change in trim angle of attack from i T = 0 to i T = -0.1 
is more than twice that from ij = + 0.1 to irp = 0 . Also, there are evi- 
dently incidence settings between 0 and -0.1 for which there will be three 
trim angles of attack, as in the preceding discussion. 

The effect of tail thickness is shown in figure 6 and the effect on 
the lift and pitching moment for zero tail height (fig. 6(a)) is seen to 
be small for the configuration chosen. However, there is a consistent 
decrease in tail lift associated with the increase of tail thickness. 

For the high-tail case (fig. 6(b)), this effect is magnified at positive 
angles of attack, with the result that the nonlinearities due to tail 
height increase with tail thickness. It should be recalled that the tail 
thickness introduced here is of a very special type (elliptic cone tail) 
and incorporates a blunt trailing edge. The present theory would predict 
no effect of thickness on the lift and only a small effect on the pitching 
moment if the tail trailing edge were sharp. On the other hand, ordinary 
slender-body theory (i.e., with no wing- tail interference) would predict 
no effect of elliptic thickness on either the lift or the pitching moment, 
regardless of whether the trailing edge is blunt or sharp (see ref. 10). 

Directional stability .- In this portion of the analysis, the empha- 
sis will be on the effectiveness of the vertical tail, and we turn to a 
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somewhat more special class of wing-tail combinations as shown in 
sketch (f) although the plan forms need not be triangular. Thus we let 



Sketch (f) 

t 2 = b/s = 0. The inverse transformation of equation ( 9 ) then gives the 
vortex positions in the a plane as 



and £1 is again given by equations (6) . Now, by direct substitution 
into equation (10) with the auxiliary formulas of equation ( 0 ) , it can 
be shown that 3 


3 The actual evaluation of the imaginary part of the square root 
indicated in equation (29) is not entirely straightforward if ambiguities 
are to be avoided. Consequently the procedure is given in the appendix. 
Taking the real part of the simpler square root has been discussed in 
Appendix B of reference 11. 
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and with these expressions one can easily write the real and imaginary 

2 

parts of the sum ^ r k a k r from which the interference lift and side 

k=i 

force can he obtained directly by means of equation (15) • It can be 
seen from equation (29) that the real part of the above sum is independ- 
ent of t x , so it follows from equation (15) that the interference lift 
is independent of the size of the vertical tail, even for (B / 0. 

The lift of the wing alone is given by equation (7) and the side 
force and lift of the tail alone are given by (see ref. 10) 


Yip - iLp = rtpU 0 


2 






(30) 


This side force and lift when added to equation (7) and to the interfer- 
ence forces calculated by means of equations (15) and (29) yield the 
total side force and lift on wing-tail combinations of arbitrary plan 
form of the type shown in sketch (f) at combined angles of attack and 
sideslip. 

Calculations have been carried out in this manner to obtain the 
variations of side force and lift coefficients, based on wing area, with 
angle of attack and sideslip for several vertical tail sizes with zero 
tail height and tail incidence (hrp = irp = 0) . These calculations were 
made only for one ratio of tail span to wing span (s 1 /s Q = l) and one 
tail length ratio (d/s 0 = 6) inasmuch as the effects of these parameters 
may be surmised from the previous discussion on longitudinal stability. 

The results are presented in figures 7 and 8 in the form of force coef- 
ficients against sideslip angle for various angles of attack. It can be 
seen from figure 7 that at the very small angles of attack, the side 
force is linear with respect to the angle of sideslip. In fact the values 
of the slopes there are given by ordinary slender-body theory for the 
tail alone. At somewhat larger angles of attack, however, the curves are 
quite nonlinear and it will be noticed that the angle of attack at which 
the curve is most nonlinear increases with the tail size. It can in fact 
be shown that the most severe nonlinearities occur at such an angle of 
attack and sideslip that one of the rolled-up vortices from the wing 
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strikes the tip of the vertical tail, since in that situation the vertical 
fin experiences a maximum side force due to the vortex. The angles of 
attack and sideslip corresponding to this condition are given, according 
to the present theory, by equations ( 6 ) if we set or £ 2 equal to 
it Q at x x = d. Thus we find 


3 = 


± 


Tt/k 


(^q/ S q) + (^t/ s q ) 

(d/s G )(l - 2 /rt 2 ) 


(31) 


In order to facilitate the fairing of the curves in figure 7 , the 
angle of sideslip given by equation ( 31 ) was calculated for the cases 
shown in the figure, and the corresponding side forces were computed. 

It will be noted that if the angle of attack is less than that of equa- 
tion ( 31 ) , then one vortex strikes the trailing edge of the vertical tail 
at the angle of sideslip indicated in equation ( 31 ), thus causing a sharp 
peak in the variation of C y with 3 . On the other hand, at angles of 
attack higher than that of equation ( 31 ), the vortex passes above the 
vertical tail and the variation of side force with angle of sideslip is 
smooth although nonlinear. In figure 7 ( 8 ) we see a very sharp peak in 
the curve for a = 0.20 which corresponds almost exactly to the angle of 
attack indicated by equation ( 31 ) for the vertical tail size t Q /s x = 1 . 0 . 

Another important point regarding figure 7 can be seen by observing 
the variation of the slope (dCY/ 83 )p_Q with angle of attack, since this 

quantity gives a measure of the directional stability; that is, since 
the vertical tail lies behind the probable center of gravity of the air- 
plane, a positive Cy^ corresponds to a negative C n p which is unstable. 

Thus it can be seen from the slopes through the origin in figure 7 that 
increasing the angle of attack has a destabilizing effect for angles of 
attack below the value given in equation (31) • It should be mentioned 
here that, according to slender -body theory (no interference), the sta- 
bility derivative (^C Y /^ 3 )p =0 is independent of the angle of attack. 

The variation of lift coefficient with angle of sideslip is given in 
figure 8 for various angles of attack. As was mentioned earlier, the lift 
is independent of the size of the vertical tail. It is observed that the 
lift is a minimum at 3=0 due to the fact that the horizontal tail is 
then in a maximum downwash field. Again, the effect of sideslip is felt 
only at the tail, so it can be concluded that at 3=0 the pitching 
moment has a positive maximum; hence sideslipping produces a nose-down 
pitching moment. 

In order to see a little more clearly the effect of angle of attack 
on directional stability, approximate calculations have been made of the 
derivative (<3c n /d3)p =0 at various angles of attack. It can be seen from 
equation ( 2 ) that a detailed calculation of the yawing-moment variation 
and subsequent differentiation with respect to 3 would be virtually 
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impossible by analytical means in view of the complicated expression of 
equation (29). Consequently, rather than performing laborious numerical 
calcinations and measuring slopes graphically, we shall simply assume 
that the center of pressure of the tail is unaffected by the wing vortices 
and lies always at the 2/3-root chord of the tail. 4 Hence, since the 
flat-plate wing experiences no side force, the derivative C no is given 


approximately (for c'/c = 1 and d/s 0 = 6) by 


' n 3 


^np - (dC n /d 0 )p_Q k - (d/2s 0 ) (dCy/dp) p_Q ~ -3(dCY/dp) p = Q ( 32 ) 

where the yawing-moment coefficient is based on the wing span 2s 0 and 
the slopes (dCy/dp)p_ 0 are measured on figure 7. The resulting varia- 
tions of C n with angle of attack are shown in figure 9 and it can be 
seen that an^increase in vertical tail size provides substantial increases 
in directional stability but that the directional stability diminishes 
considerably with increasing angle of attack. For the cases calculated 
( s i/s Q =1 ), it is observed that at about 10° angle of attack the sta- 
bility contributed by a vertical tail the same size as one panel of the 
horizontal tail (to/s! = 1.0) is equal to the stability contributed by 
a vertical tail only half that size at a = 0. 


As was mentioned earlier, a negative value of C n indicates direc- 
tional instability, but in no case has a negative value been shown in the 
above figures. It might be well to consider now the likelihood of such 
an instability. First it will be recalled that the strengths and posi- 
tions of the rolled-up vortices were calculated under the assumption that 
the circulation distribution at the wing trailing edge is elliptic, as 
predicted by slender-body theory for low angles of attack. In order to 
gain some idea of the effect of changes in circulation distribution, it 
was assumed that the circulation distribution tends toward a triangle at 
the higher angles of attack, as does the span loading at p = 0 according 
to low-speed experiments on triangular wings. Thus, keeping the same 
lift on the wing as before (eq. (7)), with a triangular circulation dis- 
tribution, the lateral spacing of the vortices becomes equal to s Q and 
the strengths then become in the same manner as equation (7), 


r = «u 0 os 0 (33) 

4 As a rough check on this approximation, the pitching moments were 
calculated in this manner for several span ratios with h>p/s 0 = 0.9568, 
d/s Q = 6 and were compared with those shown in figure 3» The average 
error for these cases was of the order of 10 percent for the angle range 
calculated. For shorter tail lengths, the error will of course be larger. 
It should be noted that with this approximation the dependence of the 
moments upon plan form is contained in the choice of the center of pres- 
sure of the tail. The choice of the 2/3-root chord implies a triangular 
tail. 
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Using this strength and spacing, then, the variations of Cy and Cl 
with p for t Q /s 1 = 1 were recalculated at a = 0.24 and the results 
are given in figures 7(b) and 8. It can be seen that this change in 
circulation distribution has magnified the effects of sideslip and has 
produced a decidedly unstable side force slope at 3=0. The effect on 
the variation of lift, however, is small. 

On an actual airplane or missile, a long fuselage whose nose is well 
ahead of the wing would further reduce the directional stability, although 
no calculations of this effect can be made here because of the complicated 
influence of the fuselage on the positions of the vortices at the tail. 
Thus it can be seen that the curves presented in figures 7> and 9 are 
probably not unduly pessimistic. 

As was mentioned previously, the greatest loss of directional sta- 
bility due to wing- tail interference is suffered when the angle of attack 
is such that at some angle of sideslip one vortex strikes the tip of the 
vertical tail. The sign of this angle of attack can be reversed simply 
by inverting the tail assembly as shown in sketch (g) . With this arrange- 
ment, the loss of directional stability can more easily be moved out of 





the flight range of angle of attack. The calculated curves of C n p 

against angle of attack are shown in figure 9 and are seen to be reflec- 
tions about the a = 0 axis of the curves obtained with the vertical tail 
above the horizontal. The arrangement shown in sketch (g) has obvious 
disadvantages in other respects, of course, as for example a negative 
dihedral effect and complications of the practical landing gear problem. 
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Lateral stability .- In this portion of the analysis we shall be 
concerned, with the variation of rolling moment with sideslip angle at 
various angles of attack. According to equation ( 3 ), the interference 
rolling moment in steady straight flight is 

1 m 

Li* = | PU 0 R j F'd(C£)- | pU Q R J F'd(SC) - pR J R ^r k a kr dx 
x=Z x=Xrpg + TE k=i 

( 34 ) 

and it is seen that the real part of the integral £ F'd(££) is to be 

evaluated at x = l and at x = x^Tg + . It must be noted here, however, 
that the contour integral is to be performed round a contour enclosing not 
only the airplane cross section but also the stream surfaces forming the 
vortex sheets. This is because of the development in references 10 and 11 
where (a) the boundary condition required that the contour of integration 
be composed of stream surfaces and (b) it was decided to enclose al 1 the 
external vortices within the contour in order to simplify the treatment of 
the analytic integrals. Here, since the nonanalytic integrand contains 
logarithmic singularities outside the body, there w ill be a contribution 
to the integral along the branch cuts representing the vortex sheets. 

To illustrate the application of equation (34), we shall confine our 
attention to flat-plate wings and symmetric cruciform tails, as shown in 
sketch (h). The plan forms, however, need not be triangular. In the 
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present case, since the external vortices arise from a separate wing, 
there is a vortex sheet common to both vortices and across which the 
potential differs by a_constant, namely the circulation T. Therefore, 
since the sign of d( ££ ) is opposite on the two sides, the contribution 
to the nonanalytic integral along this cut, regardless of its shape, is 
found to be 



vortex 

sheet 


R J - Td(y2 + z2) 

-r(y x 2 + 2 .Z - y 2 2 - z 2 2) = -r( yi 2 - y 2 2) 


(35) 


since from equations (4) and (6) = z 2 . 


As for the remaining contour integration round the airplane cross 
section, it will be convenient to perform that in the transformed a 
plane where on the circle aa = r Q 2 . The transformation for the symmetric 
cruciform tail of sketch (h) is obtained from equation (8) by setting 
t 1 = t 2 = s so that 


£ 2 = a 2 



(36) 


and r Q = s/s/2 . It follows, then, that on the circle where a = r Q e 

ig*. 

= ±4r 0 2 sin 20 d0 


ie 


(37) 


Now, recalling- that 
to the shed vortices and 
as (see sketch (b), page 


F' 


i T 

2* 


Zn(a 


°i) 


F' is the additional complex potential due 
their images , we can write F 1 in the a plane 

9) 

- Zn(a - - ln(o - cr 2 ) + ln(a - 


(38) 

and it is seen from equations (37) and (38) that the contribution of the 
starboard vortex and its image to the real part of the contour integral 
round the tail cross section is made up of integrals of the type 


R 


F'dUO s J 1 = ± 


2rr Q 2 

it 


tail 


Inf a - a. 

. \ 


Zn a 



sin 2 9 d0 

(39) 


or, since the real and imaginary parts of a on the circle are given by 
t = rocos 9 and q = rosin 9, it is found that 
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Ii = ± 


2rr 0 2 


tan 


sin 0 - — i 


.cos 0 - 


-tan 


-1 1 


sin 0 - 


cos 0 



sin 20 d0 ( 4 o) 


where - r ]_cos 0 1 and tj 1 - r^sin 0 1 are the coordinates of the star- 
board vortex in the cr plane. The entire integral round the circle is 
made up by choosing the proper sign of equation (37) and therefore of Ii 
appropriate to each half-quadrant and combining the resulting definite 
integrals. Since each half-quadrant of the circle corresponds to a partic- 
ular surface of the cruciform tail, the signs are easily chosen by observ- 
ing that d(££) = d(y 2 + z 2 ) and deciding whether this quantity is positive 
or negative in each half -quadrant as the contour of the cruciform tail is 
circumscribed in the counterclockwise sense in the physical plane. The 
resulting expression for the real part of the integral round the tail due 
to both vortices is found in this manner, after integration by parts, to be 


t</Vd(ec) = 


u j 

tail 



l<2 tan 


-1 


1 

1?! 
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(— 
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ro 4 

“3it 


A r Q 2 l 



r Q 2 \ r Q V 


Zn 


l + Ii! + 4li^ 

r o 4 


1 + 


Ii. 

r o 4 


-4 


TlT). 

ro 2 


where 



T i = R(di)i T 3. = t 2 = R (cr 2 ); r\ 2 =l(a 2 ) 


ri 2 = [ R/( Ci) ] 2 + [ I(a a )] 2 j r 2 2 = [R(a 2 )] 2 + [ I(a 2 )] 2 


and the arctangents are taken to lie between 0 and 2 it. Note that the 
expressions of equations (35) and ( 4 l) both vanish at x=x IE+ so that 

Rjp F d(^) = 0 ; also, at x = Z, yi 2 - y 2 2 = -«s 0 pd from equation (5). 

X=X TE+ 
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Turning now to the last term of equation (3*0, if we rewrite the 
transformation of equation (36) in the form 



since r Q 2 = s 2 /2, it follows directly that the vortex positions in the 
a plane are given by 

*k = |(7s k 2 - s2 + / + s2 ) ( 4 3) 


Hence the quantity a^ r appearing in the summation of equation (3*0 is 



il W 

so that the sum is given by 

^ r k a k r = r i a ir + r s a 2r = r 
k=i 

= r ( nJH - s2 + ilT^! 2 + S 2 -R^ 2 2 - S 2 - il JT7 + s 2 ) 

(45) 


Thus if we consider only zero tail incidence^ i T so that av = arj = a, 
and recall that for steady straight flight R = U Q P + iU 0 a, the last 
integral of equation (3*0 can be written in the form 



The last term of this expression is the contribution of the integral between 
the trailing edge of the wing and the apex of the tail; that is, where 
s = 0. The vortex positions £ j_ and £2 and the circulation T are given 
by equations (6) and (7), respectively. 
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The integrals appearing in equation (46) , which is valid for any 
plan for m with straight trailing edge, are of precisely the same form as 
those encountered in equation (21) for the pitching moment and may he 
treated in an analogous manner (for triangular plan forms) by making use of 
the Appendix. In the present case the second integral of equation (46) 
has a discontinuity of the real part of the square root, and the imaginary 
part vanishes, when the vortex passes through the vertical tail. Thus 
sketch (d) , page l6, is applicable to such cases if we interchange the 
labels on the two curves. The condition for the vortex passing through 
the vertical tail is yi = 


= 0 for 


by equation (6). Thus if 
gular plan forms) 


yi = 0 for 


*1 - 


< Sl 


where yi and zi are given 


hp = irp = 0, this condition becomes (for trian- 
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(47) 


If equation (47) is satisfied, then the second integral of equation (46) 
must be divided in a manner exactly analogous to that shown in equa- 
tions (23) to (26a). 


Finally, the interference rolling moment for the wing-tail combination 
of sketch (h) is obtained by substituting equations (35) , (4l), and (46) 
into equation (34). The coefficient based on the wing span 2s 0 is 

given by this result divided by pU 0 2 S w s 0 . Now, since the rolling moment 
of the symmetric cruciform tail alone is zero (see ref. 3) > and the rolling- 
moment coefficient of the isolated triangular wing is (see ref. 10) 

Cl w =-j°£ C^8) 

the total rolling-moment coefficient for the combination of sketch (h) is 
found by adding this to the above result for Cj 

Calculations have been carried out to determine the rolling moments 
for a wing- tail combination of triangular plan form having zero tail height, 
a tail-length ratio d/s Q of 6, and a span ratio s i/ s o of 1. The varia- 
tion of Ci with p was calculated at two angles of attack and the result- 
ing curves are presented in figure 10. It can be seen that at both angles 
of attack the variation of rolling moment with angle of sideslip is quite 
nonlinear and has a sharp peak at the angle of sideslip at which the star- 
board wing vortex strikes the vertical tail according to equation (31). 
However, the peaks are in opposite directions so that the curves cross one 
another. The reason for this behavior lies in the vertical location of the 
impinging vortex relative to the axis of the tail. At a, = 0.1, the vortex 
striking the tail lies close to the tail axis and therefore, due to its 
direction of rotation, induces a large negative rolling moment at the 
critical sideslip angle. On the other hand, at a = 0.2, the vortex strikes 
the vertical tail near its upper tip and therefore induces a positive roll- 
ing moment because of the velocity gradient associated with the vortex. 

For purposes of comparison, the rolling moments of the wing alone from 
equation (48) are also shown on figure 10 and it can be seen that the 
rolling moments are heavily influenced by wing-tail interference. 
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Cruciform Wing and Cruciform Tail 


The method used in the foregoing analysis of wing-tail combinations 
having plane wings could equally well be applied to those having cruciform 
wings if one had analytical expressions for the positions of the rolled-up 
vortices behind such wings. For the special case of an equal-span cruci- 
form wing in steady straight flight at 45° angle of bank, these vortex 
positions have been determined analytically in reference 13 . Consequently, 
the longitudinal characteristics of a wing-tail combination of the type 
shown in sketch (i) (the so-called " interdigitated" arrangement) can be 
studied in exactly the same way as was the arrangement shown in fig- 
ure 1 (b), if equations ( 6 ) are replaced by the expressions of reference 13 
or reference 3 for the positions of the four rolled-up vortices. 



(Note, however, that ref*. 13 uses wind axes.) It has been shown in 
rei erence 3 that the four vortices are of* ecjual strength given by 

r = */2 U 0 s oaw ( 49 ) 

Now, if we recognize once again that the vertical tail has no influ- 
ence on the purely longitudinal characteristics, a comparison of the base 
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cross section shown in sketch (j) with that of sketch (c), page 11, will 
show that equations (ll) to (l4) are directly applicable to the present 
case. Thus by simply extending the sum- 
mation of equation (12) over four rather 
than two symmetrically placed vortices, 
we find, since a 2 = ~a 1 and a 4 = -o 3 
that 


4 

l 

k=i 


r k CT k r 


-r(„ 


i r + a ir + ° 3 r + ° 3 : 




= 2TR^a lr + cT3 r ^ 


(50) 


Hence the counterpart of equation (l4) 
is found to be, after simplification 
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I 

k=i 
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r k a k r = _ , 

1 -(b/s) 


R 


f(^x + 5 S )- 







( 51 ) 


The interference lift and side force can now be obtained directly from 
equation ( 15 ) from which it follows that the interference side force 
vanishes since the above sum is real. At this point it should be noted 
that the quantity R(^ + £ 3 ) is a constant behind the wing, being equal 
to the lateral distance between the centroids of vorticity of the two 
halves of the wake. In fact, since the circulation distributions on the 
two component wings of a cruciform wing banked ^5° are identical ellipses 
(see ref. 3)> it follows that 


R (5 X + i 3 ) 


TCS 


O 


2^2 


( 52 ) 


With this relation, if we note that immediately behind the wing s = b = 0, 
we can substitute directly into equation (l6) for the interference lift. 
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Thus, making use of equation (5l) , we find 
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(53) 


and the positions 


£i z / s o and ^ 3 / s o 


can be obtained directly from 


figure 2 of reference 3 tor any given tail length. Now since the lift 
of the cruciform wing alone is the same as the lift of a plane wing of 
the same span (see ref. 3), it follows that the total lift coefficient 
based on the area of the wing, for the cruciform wing-tail combination 
of sketch (i) is, for any plan form with straight trailing edge. 





Calculations have been carried out to determine the lift curves for 
several combinations having zero tail height, zero tail incidence, and 
zero tail thickness using equation (54). Furthermore, the aspect ratio 
of the wing was taken as 2, so that the only parameters investigated here 
are the tail-span-to-wing-span ratio s 1 /s Q and the tail-length ratio 
d/s D . The resulting lift curves are presented in figure 11 for a fixed 
tail length with several values of s x /s 0 and in figure 12 for a fixed 
ratio s 1 /s Q with several tail lengths. It can be seen that the general 
character of the lift curves is very similar to that of the plane wing- 
tail combinations with tail heights different from zero. Here again 
there is a critical angle of attack at which the vortices strike the tail 
trailing edge, as can be seen by examining sketch (i), page 32. There 
are several important differences in this respect, however, between the 
characteristics of the plane and cruciform combinations . Due to the 
orientation of the vortices at the wing trailing edge with respect to 
the tail of the cruciform combination, depends on s^Sq as well 

as on h,p/s 0 and d/s 0 and no tail height is required for the existence 
of a critical angle of attack. On the other hand, the change in lift- 
curve slope at the critical angle is not so great for the cruciform since 
only half of the wing vortices strike the tail there. Also, in figure 11 
since h^i = 0, there is complete symmetry with respect to angle of attack 
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and we find a negative critical angle equal to the positive one. At the 
negative critical angle, the upper two vortices strike the horizontal 
tail whereas at the positive critical angle it is the lower two vortices • 

A final point that seems worthy of mention in this discussion con- 
cerns the existence of a so-called "leapfrog” distance for the cruciform 
wing at 45 bank. This is the distance behind the wing at which the upper 
two vortices pass between the lower two due to their mutual induced 
effects . It is clear that if the tail location and span are such that 
the vortices strike the tail trailing edge at the leapfrog distance d=dL, 
a maximum effect may be expected. According to equation (60) of refer- 
ence 3> the tail length for this condition is 


4.664 llL = 

s o 


(55) 


Hence the critical angle of attack for this maximum effect 

9.328 

a ' cr “ ~r 7 T 
jt(d/s Q ) 


is given by 


(56) 


Also from figure 2 of reference 3 it is found that if all four vortices 
are to strike the tail, the ratio of tail span to wing span s x /s 0 would 
have to be not less than 0.9 and the tail height ratio h T /s Q would be 
equal to about 1.9* It should be pointed out that for the cruciform 
wing-tail combinations for which most of the calculations were made here 
(d/s Q = 6) the value of predicted by equation (58) is probably well 

beyond the range of validity of the present theory (oe r « 28°). However, 
for the case d/s Q = 12 (see fig. 12) if the tail were raised so that 
h>p/so = 1*9> we should expect that the gentle dip in the lift curve shown 
at a = 0.06 would be replaced by an abrupt dip at a = 0.247 = 14.2°. 

Due to the complexity of the expressions given in reference 3" for 
the vortex positions and £ 3 , a direct application of equation (2) 
to obtain analytical expressions for the pitching moment appears to be 
virtually impossible. Consequently, rather than to embark on a program 
of numerical calculations, it was decided to approximate the pitching- 
moment variations by assuming that the center of pressure of the tail 
remains at its 2/3-chord position in spite of the vortices. Hence, 
since the pitching moment of the wing alone about its 2/3-chord point is 
zero, the total pitching -moment coefficient of the wing-tail combination 
about that point is given approximately by (see footnote 4-, page 25) 


* - C LT’(t + 3 ' 3 c") 

where the pitching-moment coefficient is based on the wing area and wing 
chord and the coefficient Cj^' refers to the lift on the tail in the 
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presence of the wing wake. The resulting pitching -moment variations are 
presented in figures 13 and lb and it is seen that they are comparable 
with the variations for the plane wing-tail combinations having varying 
tail span. The above discussion of the lift curves cam in fact be 
extended to the pitching-moment curves in comparing the characteristics 
of the cruciform interdigitated wing-tail combinations with those of the 
plane wing-tail arrangements with high tails. 

The influence of tail length on the value of the critical angle of 
attack and on the values of the pitching moment at that angle are clearly 
shown in figure 13 • It can be seen that doubling the tail length reduces 
the critical angle of attack by one half, but does not change the value 
of the pitching moment at the critical angle. As a matter of fact, it 
can be shown that, for given values of s 1 /s Q , b«p/s Q , and irp, the pitching 
moment of equation (57) is a function of the parameter a(d/s Q ) only. 
Hence, if we plot C m vs. a(d/s 0 )> all the points calculated for 
figure 13 will fall on a single curve; therefore, figure l4 is plotted 
in this manner. Figure l4 shows the effect of span ratio s±/s 0 and 
it is noted that the maximum effect at the critical angle is experienced 
for s±/s 0 = 0.6. However, due to the change in relative vortex posi- 
tions with distance behind the wing> the value of the critical span ratio 
will depend on the tail length for cruciform wing-tail combinations. It 
can be seen from figure 14 that if the ratio of tail span to wing span is 
sufficiently large (see = 2), the nonlinearity associated with the 

critical angle of attack is diminished substantially and the pitching- 
moment curve becomes practically linear with angle of attack. The same 
statement can of course be made in regard to the plane wing-tail combina- 
tions treated earlier since the lift of a large tail surface overshadows 
the interference effects . 


PITCHING AND PLUNGING FLIGHT 


So far the cases considered in the present report have been confined 
entirely to steady straight flight. However, the expressions of equa- 
tions (l), (2), and ( 3 ) are directly applicable to quasi- stationary 
motions (i.e., to unsteady motions whose frequencies are small compared 
with the flight velocity divided by the length of the body) . In this 
section, we shall consider longitudinal motions of this type with the 
purpose of calculating the forces and moments acting on a plane triangular 
wing-tail combination which is pitching and plunging at an angle of 
attack a v . The term pitching refers to an angular pitching velocity 
(l ^ 0); the term plunging refers to a time rate of change of angle of 
attack (a ^ 0) . 


It can be seen from equations (l) and (2) that the time derivative 
is required for the calculation of the desired forces and 
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moments, so we must now express the positions and. strengths of the wing 
vortices as functions of both x and t for the motions just described. 

In order to do this, we will first assume that the strength T of each 
rolled-up vortex at a distance x 1 behind the wing is equal to the value 
of r that existed at the trailing edge when the trailing edge occupied 
that position in space. That is to say, 


r 


- 


IE 

8t 




where ti is the time required for the wing to advance a distance x-,. 
Hence t x = x 1 /U 0 and we have 



Now rrpg is given by the jump in potential A<p at the wing trailing edge 
which depends only on the angle of attack at the trailing edge. Thus for 
the pitching and plunging case, since cp does not depend on d, 

r TE = 2 s otU 0 Ow + q(c - c x )] (59) 

where c x is the distance from the wing apex to the pivot point of the 
pitching motion. Note that this expression reduces to equation ( 7 ) if 
q = 0. Now, substituting equation ( 59 ) into equation ( 58 ) we find that 
the strength of the rolled-up wing vortex is given by 

T(x,t) = 2s 0 [U 0 o w + q(c - c x ) - dx-J ( 60 ) 

It has been assumed here that q = 0. 


The vortex position remains to be expressed as a function of 
x and t, but this can be done in a manner parallel to that for the steady 
case. Thus the slope bz x /bx is given by the velocity of a two- 
dimensional vortex pair of strength r(x,t) spaced (n/2)s 0 apart. The 
vertical position z x of both rolled-up vortices relative to the tail is 
therefore (see sketch (k)) 



Sketch (k) 
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. , P 1 r(x,t) 

•1 = a w x l - d T - iT( d - x l) - / , , ■ - 

'-'o 2n(jt/2)s 0 U ( 


dx 


= aw x 3 


-1 f T q , \1 ax, 2 1 

x ~ h T " i T( d " x i) " \ 2 aw + _Cl ^J Xl " “[J - J 


(61) 


which reduces "to equation. ( 4 -) if q — cl — 0 . Now for zero sideslip the 
lateral positions of the wing vortices are again independent of x-l and 
we have corresponding to the expressions of equation (6) 



Note that the plunging motion a, introduces a quadratic dependence of 
the vortex positions upon the distance x x as shown in sketch (k) , while 
the pitching motion q causes no additional complication over the steady 
angle-of -attack problem. It should also be pointed out for the purpose 
of taking derivatives with respect to time that the only time -dependent 
quantity in equations (60) and (62) is the angle of attack ciy, assuming 
that d = 0. 


For the purely longitudinal motions being considered here , the 

inverse transformation for a thin elliptic -cone tail is given by equa- 

2 


tion 


(ID 


and the summation 



is the same as was given in 


equation (l 4 ) . The difference in the two cases, of course, (i.e., the 
steady and the unsteady) lies in the more complicated expressions of 
equations (60) and (62) for the strengths and positions of the vortices. 
Now, according to equation (l), the interference side force and lift for 
the present case are given by 



05 

as 


and. once more the interference side force vanishes due to symmetry. The derivative 

over the tail is obtained by differentiating equation (l4) with respect to time t 
tions ( 60 ) and ( 62 ) for f and £ x and noting that 



and 


I ■ 2U ° s ° d 



(64) 


In this manner, one finds, after some manipulation, that the interference lift of the plane triangular 
wing-tail combination is, according to equation ( 63 ) 


Ll * SS b + h ■ £][*/^-g( 1 -$- H] - + i <• - '*>] - 
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(65) 


where the thickness ratio b/s is again constant over the tail. Now for the isolated wing and tail, 
the lifts due to the pitching and plunging motions are additive, and the lift coefficients, based on 
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the wing area are given according to slender-hody theory by (see, e.g., ref. 10) 
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(66) 
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Consequently, the total lift coefficient, based on the area of the wing, for the pitching and plunging 
wing-tail combination of figure l(b) is, from equations (65), (66), and (67) 
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where the expression for is given by equation ( 62 ). It can be seen that the lift is a nonlinear 

function of all three independent variables a^, a, and q. Furthermore, the lift cannot be correctly 
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calculated by adding the lifts due to each of the three motions since 
equation (68) contains products of the independent variables. That is, 
the lift due to pitching and plunging for the wing-tail combination is 
not equal to the sum of the lift due to pitching and the lift due to 
plunging . 

In reference 10 formulas were given for the stability derivatives 
of first and second order for slender wing-body combinations . In that 
analysis, as in the present one, it was found that the forces and moments 
were nonlinear with respect to some of the motions considered. However, 
in the analysis of reference 10, the forces and moments depended only on 
the first and second powers of the independent variables so that it was 
convenient to define stability derivatives with all of the independent 
variables (including the angle of attack) set to zero. In this way the 
conventional stability derivatives could be obtained by a combination 
of a small number of the derivatives of reference 10 and the latter had 
the advantage of bringing out a number of useful relationships among the 
various stability derivatives . 

In the present analysis, on the other hand, the nonlinearities of 
the forces and moments including wing- tail interference are of a more 
complicated nature . The appearance of the square root involving a quad- 
ratic function of all the independent variables in the expressions for 
the lift and pitching moment indicates that an infinite number of stabil- 
ity derivatives of the type used in reference 10 would be required to 
obtain the commonly used stability derivatives. In other words, the 
forces and moments in the present analysis contain all orders of the 
independent variables, as can be seen by expanding the square root. 
Consequently, the stability derivatives in the present analysis (includ- 
ing 0^ ) will be defined as the appropriate partial derivatives evaluated 

with all the independent variables except the angle of attack set to 
zero. Thus the stability derivatives will all be functions of the angle 
of attack and one must use the derivatives corresponding to the 

equilibrium angle of attack about which small oscillations are to be 
executed if the motion is to be calculated by means of the stability 
equations. This presumes, of course, that the equilibrium condition is 
steady straight flight at an angle of attack a w at zero sideslip. 


With the above definitions, it is easily found from equations (68) and (62) that 
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The evaluation of the quantities appearing in the expression for C L . 
is actually no more complicated than the evaluation of the integral 
already treated in finding the steady pitching moment (see eq. (21)), 
since the imaginary part can be taken after the integration. That is, 
the integral containing the imaginary part can be separated into integrals 
of the type 



and 



|fd| 

■Jx 


2 

where X = B + Cl + D| as before. Both of these integrals are reducible 
to that of equation (25) and can be handled in the same way. 


It should be noted that, although the above expressions were derived 
for triangular wing and tail, equations (69) and (70) for Cl ai.d Cl 

are actually independent of plan form and are valid for arbitrary slender 
plan forms so long as the trailing edge is normal to the x axis. The 
quantity b/s then refers to the value at the trailing edge of the tail. 


Calculations have been made to determine the first-order stability 
derivatives CL a , C L( ^, and C Ld for two tail heights and for two ratios 

of tail span to wing span. The aspect ratio of both wing and tail was 
taken as 2, the tail length ratio d/s 0 = 6, and the tail incidence and 
tail thickness were set to zero (iy = b/s = 0) . The results are presented 
in figure 15 and it is seen that the variations with angle of attack may 
be appreciable, even for the tail lying in the plane of the wing. 


The two span ratios were chosen for the calculations so that one 
falls on either side of the critical span ratio (s 1 /s Q = jc/ 4). The varia- 
tions of the stability derivatives shown in figure 15 are clearly more 
severe for the larger span ratio (sx/so = 1) since, for span ratios 
greater than the critical, the trailing vortices strike the tail trailing 
edge at the critical angle of attack. It can be seen that for s x /s 0 = 1 
and Oc r = 0.2 both CL a and Cl^ jump from a negative value to a large 
positive value as the critical angle is exceeded. For the lift-curve 
slope CL a > this was already noticeable in figures 2 and 3 showing the 

steady lift curves for the same configurations. It will also be noted 
in figure 15 that for s x / s Q = 1 the derivatives CL a and Cl^ have their 

minimum values at the critical angle of attack, even for oc r = 0 (i.e., 
wing and tail in the same plane) . This is not the case with the deriva- 
tive Cl. • The discontinuities in Cl and Cl. at the critical angle 

of attack for <xq T f 0 are evidently removed and the variations with 
angle of attack considerably reduced if the ratio of tail span to wing 
span is reduced below the critical value . 


For the case of wing and tail of equal span in the same plane, CL a 
is doubled as the angle of attack increases from 0° to 12°. On the other 
hand, if the tail span is half the wing span, the increase is only 
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about 25 percent. The variation of is about 6 percent for the 

equal-span wing-tail combination and about 5 percent for s i/s Q = 0.5 
over the same angle range. 


It has already been mentioned that in reference 10 a number of rela- 
tionships were obtained among the stability derivatives for wing-body 
combinations by evaluating the derivatives at a = 0 . Although, as 
already discussed, there seems to be little point in defining stability 
derivatives for wing-tail combinations in the same manner, it is never- 
theless of interest to see whether a similar type of relationship can be 
found at =0. In fact if we set =0, it becomes apparent that 

equations (69) and (70) are then related according to 


or 
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( 72 ) 


so that Cr at ciy = 0 for the wing-tail combination can be calculated 
from the lift-curve slope of the combination at aw = 0 and the lift-curve 
s ^-°P e of the isolated tail. If there is no tail = 0 and equa- 

tion (72) reduces to the relationship given in reference 10 between Cl 
for wing-body combinations, with the chord c as the reference 

length . 


LIMITATIONS OF THE THEORY 


It has been assumed in the foregoing analysis that the vortex sheet 
leaves the wing trailing edge as a flat sheet and that it becomes fully 
rolled up ahead of the tail. The former assumption does not admit sepa- 
ration of the flow from the wing surface and therefore implies smal i 
angles of attack. The latter assumption, on the other hand, seems to 
imply a high angle of attack. Actually, it requires that the tail length 
be greater than the distance in which the wing vortex sheet becomes fully 
rolled up, that is (see, e.g., ref. 3 or 12), 



( 73 ) 
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•where k, depends on the circulation distribution. According to Kaden 
(ref. lU) the value k x = 0.28 gives the distance required for the vortex 
sheet behind a plane wing with elliptic circulation distribution to 
become "essentially rolled up." It can be seen from equation (73) that 
the restriction imposed on the configurations which can be treated becomes 
more severe as the angle of attack approaches zero. Hence the present 
calculations might be considered as a "moderate angle -of -attack" analysis 
for slender wing-tail combinations . 

In view of equation (73), it is interesting to consider other analy- 
ses in which wing-tail combinations have been treated under the assumption 
of a flat vortex sheet lying in the plane of the wing and tail. The 
requirement for such results to be applicable is clearly 



ki 




( 7*0 


The "much greater" sign is introduced because the assumption of such 
analyses is not simply that the vortex sheet is not fully rolled up at 
the tail but rather that it remains completely flat - hence the more 
stringent requirement . It is evident then that if a wing-tail combina- 
tion with a low-aspect-ratio wing is treated under the above assumption, 
the results are applicable only for vanishingly low lift coefficients or 
for extremely short tail lengths . There have recently been some such 
analyses (e.g., refs. 15 and l6) in which wing-tail combinations are 
treated with no apparent restrictions on aspect ratio or tail length. 

The results of such investigations when applied to low-aspect-ratio wings 
or sizable tail lengths should be viewed with caution. It should also 
be mentioned that the analyses of references 15 and l6 can lead only to 
stability derivatives which are independent of the angle of attack. 
Furthermore, the vortex sheet is assumed to lie in the plane of the tail 
at all (small) angles of attack. Although the vortex sheet does indeed 
lie in the plane of the tail at a = 0 (for h T = 0) , the derivatives 
calculated at a = 0 may nevertheless be in error since no account has 
been taken of the change of position of the vortex sheet with angle of 
attack. 

In reference 15 the conclusion is drawn that the lift induced on the 
tail by an angle of attack of the wing in forward flight is equal to the 
lift induced on the wing by an angle of attack of the tail in reverse 
flight. That this conclusion is not in agreement with the present analy- 
sis can be seen by examining equation ( 17 ) for the interference lift in 
steady straight flight. For the slender case, as is treated in the 
present paper, reversal of the flight direction has no influence on the 
given expression other than to interchange s x and s Q . It is clear then 
that the above conclusion agrees with the present analysis only for the 
obvious case of s x /s 0 = 1. The difference, of course, stems from the 
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different assumptions regarding the vortex sheet shed by the wing. The 
appropriateness of the different assumptions can only be checked by 
experiment . 

A further assumption of the present analysis is that the positions 
of the vortices shed from the wing are not influenced by the presence of 
the tail. Thus the vortex lines coincide not with the actual streamlines 
but with the streamlines in the absence of the tail. Hence, while the 
condition of no flow through the tail has been satisfied in the analysis, 
the vortices have nevertheless been permitted to penetrate the tail 
surface. In fact the most dramatic effects on stability are predicted 
for conditions under which the wing vortices either touch or pass through 
the tail. Although the present assumption that the tail does not influence 
the vortex positions is commonly made in calculations of wing-tail inter- 
ference, the implications mentioned above seem to warrant investigation of 
this point. Therefore an experiment was conducted to investigate the 
influence of the tail on the vortices shed by th*e wing. This experiment 
will be discussed in the next section. 


The rolled-up vortices were assumed in the analysis to be idealized 
line vortices having no viscous cores. This assumption will have no 
effect on the calculated results so long as the cores do not touch the 
tail. However, when the distance from the tail to the center of the 
vortex is smaller than the radius of the core, the viscous core will 
change the downwash distribution at the tail location and affect the 
calculated lift on the tail. In order to get some idea of the order of 
magnitude of this effect, a calculation was made by the method of reverse 


flow (see ref. 7) to determine the 
lift on a plane tail in the presence 
of two viscous vortices whose centers 
intersect the trailing edge of the 
tail as shown in sketch ( l ) . For 
this calculation, it was assumed that 
the cores rotate as solid bodies, 
giving the downwash distribution 
shown in the absence of the tail. It 
was further assumed that the cores lay 
entirely within the span of the tail 
as shown, and that the span loading 
of the tail in reverse flow is ellip- 
tic. With these assumptions, it was 
found that if the vortex centers are 
located at y = ±s 1 /2 and the core 
radius is s-l/ 4 (as shown), then the 
negative lift on the tail due to the 
vortices is reduced by 12 percent 
from that calculated with no cores . 
Now, since the centers of the vortices 
were assumed to intersect the tail 
trailing edge, the above calculation 
corresponds to a = ac r . Furthermore, 
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since the viscous cores can have no effect on the tail lift if they do 
not touch the tail, it can be concluded that the over-all effect of the 
viscous cores on the calculated lift and pitching-moment curves of fig- 
ures 2 to 6 will be simply to round off the sharp peaks at a = acr- 

Although the present analysis has been restricted to slender wing- 
tail combinations for which the vortex sheet leaves the wing as a flat 
sheet, it is important to note that the origin of the trailing vortices 
is immaterial to the method of calculation of wing-tail interference. 

If, for example, the wing is a delta wing with sharp leading edges, there 
may be separation vortices above the wing. These can be handled by the 
present method, provided that their positions and strengths over the wing 
are known. Similarly, if the wing is a sweptback wing (not necessarily 
slender), the vortices from the wing may be caused to roll up quickly 
because of tip-stalling or other viscous phenomena. Nevertheless, if 
the tail is slender and the vortices are rolled up ahead of the tail, the 
techniques used in the present analysis are still applicable for the 
calculation of the interference forces and moments. As a matter of fact, 
pitching -moment curves having the same general nonlinear character as 
those of figures 2 and 3 have been measured on swept-wing airplane models 
having high tails (see, e.g., refs. 17 and l8) . Some of the trends in 
directional stability which have been calculated in the present report 
are also evident in the data of reference l8. 


WATER-TANK EXPERIMENT 


In order to investigate the influence of the tail on the paths of 
the vortices shed by the wing, a simple model was constructed consisting 
of an identical triangular wing and tail mounted on a thin plate, as 
shown in sketch (m) . The plate acted as a model support and was attached 
to a rack which was driven vertically into a water tank by a gear mounted 
on an electric motor. The wake of the wing was made visible on the sur- 
face of the water by means of white poster paint applied to the beveled 
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trailing edge of the wing. As the wing passed through the water surface, 
the paint remained on the surface and the resulting wake patterns were 
photographed from above by a motion picture camera. The wing and tail 
were flat plates with rounded leading edges and beveled trailing edges, 
and the central plate was cut down to minimize its influence on the wing 
wake (see sketch (m)). The tail was mounted at zero incidence (iip = 0) 
in a tail-high position (hrji/s 0 = 0 . 96 ) so that at a critical angle of 
attack (ac r = 0 . 2 ) the analysis will predict that the wing vortices inter- 
sect the tail at its trailing edge. 

The model was driven into the tank at three angles of attack, one 
below, one above, and one at the calculated critical angle of attack, 
and selected frames of the motion pictures obtained are presented in 
figure l 6 . The vertical and horizontal lines on either side of the frame 
are reference markers which are fixed relative to the tank (i.e., relative 
to the flight direction) . The first picture in each case shows the vor- 
tex sheet leaving the wing trailing edge as a flat sheet . 5 At this point , 
the wing trailing edge has just penetrated the water surface, and the 
strut connecting the wing and tail can be seen above the water surface. 

For a = 12° and 16° the tail is also visible but is out of focus , being 
closer to the camera. The second picture in each series shows the tail 
entering the water surface, and the remaining pictures show the develop- 
ment of the wing wake at various positions over the tail until the last 
frame shows the wing wake at the trailing edge of the tail. 

In figure 16(a) there seem to be no obvious effects attributable to 
the tail, but the mounting strut has a marked effect on the shape of the 
center of the vortex sheet. It can be seen that the upward sweep of the 
center of the sheet follows the retreating edge of the strut. 

Figure 16(b) shows the development of the wing wake over the tail at 
the calculated critical angle of attack. According to the analysis in 
the foregoing sections of this report, the wing vortices should intersect 
the tail at its trailing edge at this angle of attack. It was in fact 
largely this prediction that prompted the present experiment, in view of 
the discussion of the previous section. It is apparent from figure 16(b) 
that, although the vortices are quite distorted at the trailing edge of 
the tail, the centers of the vortices do very nearly coincide with the 
tail trailing edge. The tail seems to split the viscous cores of the 
vortices at the critical angle of attack. The asymmetry of the center 
portion of the vortex sheet seen in the photographs is probably due to 
an asymmetric disturbance produced by the strut. 

In figure 16(c) a higher angle of attack is shown and the tail appears 
to sever the vortex sheet and produce considerable distortions above the 
surface of the tail. 

Subsequent detailed investigation has shown the existence of a rela- 
tively weak pair of separation vortices above the wing surface at all of 
the angles of attack studied here. 
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The distances behind the wing as shown on the above-mentioned pic- 
tures were measured by means of a tape which moved with the model and 
recorded on the film the distance between the wing trailing edge and the 
water surface. Thus it was a simple matter to obtain quantitative measure- 
ments of the vortex paths , provided the centroids of vorticity of the 
vortex sheet could be defined. Actually it is not possible to define these 
positions accurately from photographs, but it was assumed that the lateral 
position of the centroid of vorticity of each half of the wake could be 
taken as halfway between the outermost winds of the rolled-up portion of 
the sheet. The results of such measurements are shown in figure 17 as 
the variation of vortex span with distance behind the wing. The position 
occupied by the tail is indicated in the figure and it can be seen that 
the influence of the tail at all three angles of attack is such as to 
decrease the vortex span at the trailing edge of the tail. 

The theoretical asymptotic vortex span of (jr/4)s 0 for the wing alone 
is in close agreement with the experimentally observed positions at the 
lower angles of attack until the influence of the tail is felt at the 
station containing the apex of the tail. The sizable inward shift ahead 
of the tail at the largest angle of attack is associated with the change 
in circulation distribution on the wing as the angle of attack is 
increased. Inasmuch as the interference lift calculated in the analysis 
depends on the positions of the wing vortices at the trailing edge of the 
tail, it would be of interest to calculate the effect on the interference 
lift of the observed inward movement due to the tail. However, the inter- 
ference lift is a function of the vertical positions as well as the lateral 
positions of the vortices, and it was felt that the vertical positions of 
the centroids of vorticity could not be defined with any degree of accuracy 
because of the unrolled-up portion of the vortex sheet. It is interesting 
to note that at the critical angle of attack Zj_ = 0 and the interference 
lift does not depend on the lateral spacing of the vortices so long as 

they lie within the span of the tail. (This point has been discussed in 

the section titled "Longitudinal Stability.") Consequently, at the criti- 
cal angle of attack the influence of the tail on the vortex positions has 
no effect on the interference lift for the span ratio tested. For other 
angles of attack or for sufficiently smaller tail spans this will not be 
the case . 

The systematic tendency of the tail to draw the vortices together as 
observed in figure 17 suggested the possibility that the water-tank obser- 
vations have been influenced by surface tension effects. Therefore, 
observations were made of the three-dimensional wake patterns below the 
surface of the water to see whether any significant changes were occurring 

at the water surface due to surface tension. This was done by thinning 

the poster paint so that some of it would adhere to the wing after the 
latter was submerged, thus extending the visible wake patterns below the 
water surface . Observations were made at various angles of attack, and 
it was found that no surface tension effects were discernible as long as 
the angle of attack was greater than 6 . Below that angle of attack, 
the difference between the patterns on the surface and below it became 


MCA TN 3725 


51 


visible . Apparently , the criterion of a low Weber number (ratio of sur- 
face tension forces to dynamic forces) was satisfied for angles of attack 
above 6°. 

The assumption made in the analysis that the vortex sheet becomes 
fully rolled up ahead of the tail appears to be substantiated in the case 
of the model tested here. 


CONCLUDING REMARKS 


A theoretical investigation has been made of some of the aerodynamic 
characteristics and stability problems associated with slender wing-tail 
combinations. The vortex sheets shed by the wing have been assumed to be 
fully rolled up at the tail and to follow the same paths as calculated in 
the absence of the tail. 

Mathematical expressions have been derived for the interference 
forces and moments acting on the tail. From these, calculations have been 
made of the effect of changes in tail height, tail length, tail incidence, 
tail thickness, and ratio of tail span to wing span. The calculated 
variations of forces and moments with angles of attack and sideslip were 
found in some cases to be highly nonlinear. Changes in the height of the 
horizontal tail and in the ratio of tail span to wing span were found to 
have a pronounced effect on static stability as well as on the stability 
derivative Cl - . The calculated results indicated a definite critical 
angle of attack for a given wing-tail combination and a critical ratio of 
tail span to wing span at which abrupt changes in the aerodynamic charac- 
teristics occurred in certain cases. In general, the most drastic effects 
are predicted when the vortices shed from the wing strike or pass near the 
tips of the tail trailing edge. 

A water-tank experiment was conducted in order to observe the 
behavior of the wing wake in the presence of the tail. A plane -wing-tail 
combination with a high tail was tested and photographs indicated that 
the tail caused the wing vortices to be shifted inboard appreciably for 
the case tested (wing and tail of equal span). The vertical positions, 
however, were apparently relatively unaffected, and the tail was observed 
to sever the vortex cores at the calculated critical angle of attack 
(within the experimental accuracy) . 

The influence on the calculated results of the assumptions made 
regarding the vortex wake have been considered, and a conclusion drawn 
on the basis of a flat vortex sheet in the plane of the tail was compared 
with the results of the present analysis. The choice of such assumptions 
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must in the last analysis be made on the basis of experiment, and for the 
cases treated in the present paper the assumption of the fully rolled-up 
vortex sheets seems justified, provided that the tail length is not 
extremely short. 


Ames Aeronautical Laboratory 

National Advisory Committee 
Moffett Field, Calif., 


for Aeronautics 
May 21, 1956 
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APPENDIX 


EVALUATION OF THE IMAGINARY PART OF A COMPLEX SQUARE ROOT 


A complex square root of the 
and c is real has branch points 
factored form 

7 = -id - ic)(5 + ic) 


Now each factor can be written in 
polar coordinates referred to one of 
the branch points. Thus (see 
sketch (n) ) let 

, i<p f . • 

c - 1C = P x e ; c, + 1 C = p 2 e 

where cp and V are both limited to a 
range of, say, -(«/2) to +(3«/2) (to 
give the proper sign changes through 
the line segment shown) . This 
enables one to write 


type + c 2 where £ = y + iz 

£ = ±ic and can be written in its 


z 



/cp+V\ 

J? * +1 sln(2|i)] 


so that the imaginary part is 


where 




P x = |*)y 2 + (z - c ) 2 

p 2 = |7y 2 + ( 2 + c) 2 
cp = tan -1 
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t = tan -1 



The sign of the imaginary part of the square root is therefore determined 
by the quadrants of cp and \|r which depend on the position of the point £ 
relative to the branch points £ = ±ic and also on the range of 9 and \|r 
specified above. 
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Figure 1.- General orientation of wing and tail. 
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(b) 


Type of configuration treated for longitudinal stability analysis. 


Figure 1.- Concluded. 
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Figure 2.- Effect of horizontal tail height on the variations of lift 
and pitching moment with angle of attack; plane triangular wing 
and tail of aspect ratio 2, d/s 0 = 6, sJsq = 1, i t = b/s = 0. 
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Figure 3.- Effect of tail-span-to-wing-span ratio on the variations of 
lift and pitching moment with angle of attack; plane triangular wing 
and tail of aspect ratio 2, d/s Q = 6, h^/s Q = 0 . 9568 , i^ = b/s Q = 0. 
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(a) h t /s Q = 0 

Figure Effect of tail length on the variations of lift and pitching 
moment with angle of attack; plane triangular wing and tail of aspect 
ratio 2, Sl /s 0 = 1, i t = b/s Q = 0. 


62 


MCA TN 3725 




(b) h t /s 0 = 0.9568 
Figure 4.- Concluded. 
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Figure 5«- Effect of tail incidence on the variations of lift and 
pitching moment with angle of attack; plane triangular wing and 
tail of aspect ratio 2 , Sj/sq = 1 , b/s 0 = 0 , d/s Q = 6. 


64 


NACA TN 3725 




NACA TN 3725 


65 




Figure 6.- Effect of tail thickness on the variations of lift and 
pitching moment with angle of attack; plane triangular wing and 
tail of aspect ratio 2, Si/s 0 = 1> it = d/s 0 = 6. 
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(b) h t /s Q = 0.9568 


Figure 6.- Concluded 
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(a) t Q /s x = 0.5 

Figure 7.- Effect of angle of attack on the variation of side force with 
angle of sideslip; plane wing with horizontal and vertical tail, 
s i/ s o = 1 > d / s o = 6, h t /s Q = i t = 0. 
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(b) t 0 /s x — 1.0 


Figure 7»- 


Continued 
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(c) to/s-L = 1.5 
Figure 7 • - Concluded . 
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Figure 8.- Effect of angle of : attack on the variation of lift coefficient 
with angle of sideslip for all vertical tail sizes; plane wing with 
horizontal and vertical tail, s^/sq = 1, d/s Q = 6 , h-^/s^^ ij- = 0. 
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Figure 9«~ Effect of vertical-tail size on the variation of yawing moment 
due to sideslip with angle of attack; plane triangular wing with tri- 
angular horizontal and vertical tail., s-^/sq = 1, d/s 0 = 6, 

■^t/ s o “ i-t = j A w = A-^ = 2. 
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Figure 10.- Effect of angle of attack on the variation of rolling moment 
with angle of sideslip; plane triangular wing and cruciform in-line 
tail, s i/ s o = > ^/ s o = ^ > -^t = ^w = ^t = 2. 
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Figure 11.- Effect of tail-span-to-wing-span ratio on the variation of 
lift with angle of attack; cruciform interdigitated wing-tail 
combination at 45° bank, d/s Q = 6, h-t/s 0 - it = 0, A w = A-fc = 2. 
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Figure 12.- Effect of tail length on the variation of lift with angle 
of attack; cruciform interdigitated wing-tail combination at 45° 
bank, bJb q = 1 , h t /s Q = i t = 0, A w = A t = 2. 
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Figure 14.- Effect of tail-span-to-wing-span ratio on the variation of pitching moment with tail 
length parameter ad/s 0 ; cruciform interdigitated wing-tail combination at 45 bank, 

h t/ s o = H = °> A w = A t = 2 * 
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Figure 15 •“ Effects of horizontal tail height and ratio of tail span to 
wing span on the variations of the lift derivatives with angle of 
attack; plane triangular wing and tail of aspect ratio 2, d/s 0 = 6 , 
i t = b/s = 0. 
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x i/s Q = 0 


x i/ s o = 4.28 




x x /s 0 = 4.84 


x i/ s o = 5.12 




x i/ s o = 5.46 


A-20904 

x i/ s o = 5.96 


(a) a = 8° 


Figure 1 6 .- Photographs of the wing wake in the presence of the tail; 
identical triangular wing and tail of aspect ratio 2 , d/s Q = 6, 
^t/ s o = 0 * 96 >j i-^ — 0. 













Figure 1 6 .- Continued. 
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Figure 16.- Concluded. 





NACA - Langley Field, Va. 



Figure 17 • - Water-tank measurements of the influence of the tail on the lateral positions of the 
vortices shed by the wing; identical triangular wing and tail of aspect ratio 2. d/s n = 6, 
V s o = 0.96, i t = 0. 
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